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^ ; ABSTRACT 

The known /? Cephei star HD 180642 was observed by the CoRoT satellite in 2007. From the very high-precision light curve, its 
pulsation frequency spectrum could be derived for the first time (Degroote and collaborators). In this paper, we obtain additional 
^ ' constraints for forthcoming asteroseismic modeling of the target. Our results are based on both extensive ground-based multicolour 

, photometry and high-resolution spectroscopy. We determine T^g = 24 500+1000 K and log^ = 3.45+0.15 dex from spectroscopy. 

OA ■ The derived chemical abundances are consistent with those for B stars in the solar neighbourhood, except for a mild nitrogen excess. 

A metallicity Z = 0.0099+0.0016 is obtained. Three modes are detected in photometry. The degree C is unambiguously identified for 
CO ' two of them: C = Q and ^ = 3 for the frequencies 5.48694 d"' and 0.30818 d"' , respectively. The radial mode is non-linear and highly 

, dominant with an amplitude in the U-filter about 15 times larger than the strongest of the other modes. For the third frequency of 

^— «j . 7.36673 d"' found in photometry, two possibilities remain: ^ = or 3. In the radial velocities, the dominant radial mode presents a 

' so-called stillstand but no clear evidence of the existence of shocks is observed. Four low-amplitude modes are found in spectroscopy 

' and one of them, with frequency 8.4079 d"', is identified as {£,m) = (3,2). Based on this mode identification, we finally deduce an 

equatorial rotational velocity of 38+15 km s"'. 
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1 . Introduction [T993I Gautschy & Saio [T993] l. Given that mainly low-degree 

rr,, „ , ■ , ^ •„ ■ low-order pressure and gravity modes are observed, these stars 

The p Cephei stars are a hoinogeneous group of oscillating BO- ^^^^^^^^ ^ ^^^.^ asteroseismic studies of 

B3 stars that have been stud^as a class for more than a cen- ^^^^^^^^^ ^^^^.^^ ^^^^^ 
tury. Stankov & Handler (I2005D compiled an overview of the 

observational properties of this group of stai-s. The oscillations I" principle, these studies would allow us to calibrate the 

of Cephei stars are explained in terms of the k mechanism that evolutionary models across the entire main sequence, although 

operates on iron-peak elements (Dziembowski & Pamyatnykh Cephei stars seem to be observed while in the second 

part of their main-sequence lifetime (Stankov & Handler 20051). 

* Based on data gathered with the 1.2m Mercator telescope Moreover, theory predicts the occurrence of stellar oscillations 

Roque de los Muchachos, La Palma, the 90cm telescope at Sierra in B stars after the main sequence (Saio et al. l2006l) and oscilla- 

Nevada Observatory, Spain, the 1.5 m telescope at San Pedro tions in such stars have indeed been found (Lefever et al. l2007l l. 

Martir Observatory, Mexico, the Im RCC and 50cm telescope at - -.j-u u i-ujfc 

, , , „ ■ o ■ r , , ^, TT Asteroseismic studies have been accomplished for 5 

the Piszkesteto Mountain Station ot Konkoly Observatory, Hungary, „^ , • , , . i-„„„ i^ ■ ^ rr., ,• 

the 2.2m ESO telescope (ESO Programme 077.D-031 1; ESO Large /^Cephei stars only (see Aerts|2008]for a review). These studies 

Programme 178.D-0361) at La Silla, Chile, the 1.93m and 1.52m tele- were based on intensive ground-based observing campaigns and 

scopes at the Haute-Provence Observatory, France led, for the first time, to constraints on the core overshoot and in- 

** Current address: Laboratoire AIM, CEA/DSM-CNRS-Universite ternal rotation of massive B-type stars. Despite these successes, 

Paris Diderot; CEA, IRFU, SAp, centre de Saclay, F-91191, Gif-sur- we clearly need to increase the number of asteroseismically de- 

Yvette, France termined models for these types of stars. This will become es- 
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Table 1. Logbook of the photometric observations of HD 180642. 



Observatory 


N 


AT 


HJD begin 


HJD end 


filter(s) 


/i 


La Silla 


20 


186.7 


50562.8 


50749.6 


UBiBBjViVG 


5.4870(2) 


La Palma 


171 


1183.8 


52421.6 


53605.4 


UB1BB2V1VG 


5.48693(1) 


SPMO 


113 


5.3 


53192.7 


53198.0 


u, V, b, y 


5.494(3) 


SNO 


222 


359.2 


53205.4 


53564.6 


u, V, b, y 


5.48693(3) 


KO 


84 


319.2 


53260.3 


53579.5 


Johnson V 


5.49(1) 


KOCCD 


261 


26.1 


53921.4 


53947.5 


Johnson V 


5.4858(3) 


ASAS 


363 


2797.6 


51979.9 


54777.5 


Johnson V 


5.486915(7) 


Stromgren 


335 


371.9 


53192.7 


53564.6 


u, V, b, y 


5.48693(2) 


U 


526 


3042.6 


50562.8 


53605.4 


U, M 


5.48695(2) 


B 


526 


3042.6 


50562.8 


53605.4 


B, V 


5.486947(2) 


V 


526 


3042.6 


50562.8 


53605.4 




5.486944(2) 


V+ 


1234 


4214.7 


50562.8 


54777.5 


V,y, V 


5.486938(5) 



The number N of datapoints retained for the analysis, the total timespan AT (in days) and the available filters are given. The last column gives the value of the 
dominant frequency (see Sect. 12. 2t expressed in d"' with its error between brackets. The HJD is given with respect to HJDq = 2400000. Stromgren refers to merged 
SPMO and SNO Stromgren data. The merged Geneva U and Stromgren u, B and i', and V and y data, are denoted as U, B, V data, respectively. The V+ dataset 
corresponds to the V data to which the KO Johnson V data and ASAS Johnson V data are also added. 



pecially important in the near future to fully exploit the data 
becoming available from ongoing dedicated ground-based ob- 
serving campaigns (see Uytterhoeven 2009jfor a review) as well 
as ongoing space missions such as MOST, CoRoT and Kepler 

The Bl.Sn-III star HD 180642 (V1449 Aql, HIP94793, 
V = 8.29 mag) was discovered as a candidate j8 Cephei star 
during the Hipparcos mission by Waelkens et al. ( I1998I ). 
Afterwards, Aerts (|2000l l used Geneva photometry to identify 
its dominant mode as a radial one, confirming it to be a/? Cephei 
star. Moreover, HD 1 80642 is among known /3 Cephei stars of the 
highest amplitude. Its peak-to-peak V amplitude in the Geneva 
filter is 0.078 mag. For comparison, the mean value for the 
/3 Cephei stars in the catalog of Stankov & Handler (2005) is 
~0.025 mag, for visual filters. 

Our studied target was observed by the CoRoT satellite dur- 
ing a long run (156 days) between May and October 2007. We 
note that it is the only known /3 Cephei star in CoRoT's core pro- 
gramme. A detailed modeling and interpretation of the CoRoT 
data is presented in Degroote et al. (2009). These latter authors 
modelled this very high-precision light curve with 1 1 indepen- 
dent frequencies and 22 low-order sum and difference combi- 
nation frequencies. Nine of the independent frequencies are in 
the range expected for /3 Cephei stars, between 5 and 9d 
Moreover, 5 of the 33 frequencies are in the regime of high-order 
g modes with frequencies below 2d"'. 

To complement the space white-light data, we conducted 
a project to collect ground-based multicolour photometry and 
high-resolution spectroscopy for selected primary and secondary 
CoRoT targets, and we refer to Poretti et al. (2007 ) and 
Uytterhoeven et al. (120081) for a description of this campaign. In 
this paper, we present the results of the analysis of the ground- 
based datasets for HD 180642. We aim to determine the funda- 
mental parameters and chemical abundances of the object and 
to derive the wavenumbers {{, m) of its low-amplitude modes, 
which are also detected in our ground-based measurements. 

The paper is divided into two parts. Sections 2 and 3 are 
devoted to photometry and spectroscopy, respectively. In both 
sections, we first describe our observations and data reduction. 
Afterwards, we describe our frequency analysis followed by our 
mode identification. In the spectroscopic part, we also present 
an abundance analysis. We end the paper with a summary of our 
results in Sect. 4. 



2. Photometry 

2.1. Data description 

The star was monitored with the 1.2m Mercator telescope 
(MER) at La Palma in the framework of a long-term photometric 
monitoring programme dedicated to pulsating stars. For our tar- 
get, we completed observing campaings of Geneva seven-colour 
(UB1BB2V1VG) high-precision photometry during two seasons 
(2002, 2005) with 80-day and 57-day durations, respectively. 
The star was observed a few times every clear night. In total, 
we used 171 7-colour Geneva measurements with a time span 
of 1184 days. Aerts (120001) had aheady obtained 20 datapoints 
of the star with the same instrument attached to the 0.7m Swiss 
telescope at La Silla. 

With the goal of gathering additional measurements with 
a higher sampling per night, HD 180642 was observed with 
the twin Danish six-channel uvby/3 photometers attached to the 
90cm telescope at Sierra Nevada Observatory (SNO), Spain, and 
at the 1.5 m telescope at San Pedro Martir Observatory (SPMO), 
Mexico. All the data were collected in the four Stromgren uvby 
filters. At SNO, a total of 222 high quality datapoints were col- 
lected during 1 night in July 2004, 5 nights in June 2005, and 6 
clear nights in July 2005. At SPMO, our target was observed dur- 
ing 4 nights in July 2004, resulting in a total of 1 13 data points. 

Johnson V measurements were also obtained with the 
50 cm telescope at the Piszkesteto Mountain Station of Konkoly 
Observatory (KO) during 2 nights in September 2004 and 3 
nights in July 2005. Since the V data showed many scattered 
points, we removed all data deviating by more than 0.05 mag 
from the fit with the dominant frequency (see Sect. |2.2| i. We re- 
tained a total of 84 datapoints for the analysis. The comparison 
star used at all sites was HD 181414 (V = 7.068 mag, A2). 

During a total of 5 nights in July and August 2006, 
HD 180642 was also observed in Johnson V filter with a 
Princeton Instruments Vers Array 1300B CCD detector attached 
to the Im RCC telescope (f/13.5) at KO. The data were binned 
such that the total integration time of the binned datapoints was 
about 3 minutes. In total, we retained 261 datapoints. Since the 
same comparison stars were monitored not every night in the 
CCD frame, we reconstructed relative magnitude differences us- 
ing common bright targets. 

Finally, we used the Johnson V ASAS publically available 
photometry from Pigulski & Pojmaiiski (i2008) . We do not in- 
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Fig. 1. Phase diagrams for the ultraviolet, blue, and visual merged datasets U, B, V+ (dots) of HD 180642 folded according to the 
dominant frequency /; = 5.48694d The full gray lines are leastsquares fits including /i, 2/i and 3/i. 
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Fig. 2. SCARGLEperiodograms calculated from the combined 
Geneva U and Stromgren u data. The uppermost panel shows 
the spectral window of the data. All subsequent panels show the 
periodograms at different stages of prewhitening. From top to 
bottom: periodogram of the observed light variations, of the data 
prewhitened with f\ and its two harmonics 2f\ and 3/i , of the 
data subsequently prewhitened with f^, and of the data subse- 
quently prewhitened with f^. The three frequencies found are 
indicated by dashed gray lines. The light gray lines indicate the 
4 signal-to-noise (S/N) level. It is computed as the average am- 
plitude over a frequency interval with a width of 2 d ' in an over- 
sampled SC A RGLE periodogram obtained after final prewhiten- 
ing. The highest peak in the third panel corresponds to a one-day 
alias (I+/2O of the true frequency (see text). 



elude the Hipparcos photomety in our analysis because the fre- 
quency value of the dominant mode (5.48709(3) d"') in these 
earlier observations, differs significantly from the values derived 
in the other photometric measurements (see Sect. 12. 2t . A sum- 
mary of the datasets used is listed in Table[T] 

2.2. Frequency analysis 

We first searched for frequencies in all filters of the individ- 
ual datasets, using the SCARGLE (Scargle [19821 ) method and 
the leastsquares power spectrum method (Vamcek |1971| ). These 



methods produced identical results within their respective un- 
certainties, so we list only the outcome obtained for the SCAR- 
GLEperiodograms. In all datasets and filters, we confirm the 
dominant frequency. The value of f\ found in each separate 
dataset is given in the last column of Table[T] where the un- 
certainty in the last digit is given in parenthese. This frequency 
uncertainty was calculated to be cry - V(6)o"std/^ V(A^)A/Ar 
(Montgomery & O' Donoghue l 1 999l , where 0-^^^^ is the standard 
deviation of the final residuals. A/ the amplitude of the fre- 
quency /, and tsJ the total timespan of the observations. 

To increase the frequency precision, we combined datasets 
for given filters whenever possible. We merged the SPMO and 
SNO Stromgren data (termed Stromgren hereafter), and con- 
structed ultraviolet, blue, and visual ground-based light curves 
(U,B,V hereafter) by merging the Geneva U and Stromgren u, B 
and v, and V and y data, respectively. To combine the data, we 
calculated the preliminary frequency solution of the individual 
subsets, checked that the amplitudes of the dominant mode in the 
Geneva/Stromgren U/m, B/v and V/y are the same within the er- 
rors, and realigned the subsets at the same mean brightness level. 
For a description of the procedure, we refer to Uytterhoeven et 
al. (2008, their Sect. 2.1). We note that in the SNO dataset it was 
also necessary to align the data obtained in 2004 with the data 
taken in 2005. 

We also constructed a more extensive visual band light curve 
(termed Vh- hereafter) by also adding the KO Johnson V data 
and the ASAS Johnson V data to the Geneva V and Stromgren 
y data. The results for the dominant frequency of these merged 
sets are also listed in Table[T] The refined value of the dominant 
frequency for each of U, B, V, and Vh- is /i = 5.48694 d"'. 

A harmonic fit to the merged U, B, V light curves is shown 
in Fig. [T] and infers the prominent presence of harmonics of fx. 
The harmonics of f\ were also derived by a frequency analysis 
of the residual light curves after prewhitening with fx . We found 
that 4/1 no longer has a significant amplitude. 

In the periodograms of the residuals of the individual and 
merged datasets, the highest power occurs at low frequency be- 
tween 5 and 9 d In the individual datasets, none of the candi- 
date frequencies reaches a high enough amplitude to be signifi- 
cant, where we adopt the criterion derived empirically by Breger 
et al. ( 119931 ) and studied further by Kuschnig et al. (199 7 ), who 
studied the variability of HST guide-stars. They concluded that a 
peak reaching at least 4 times the noise level in the periodograms 
has only a 0. 1 % chance of having been produced by noise. These 
results were confirmed by De Cat & Cuypers (120031 ) for the case 
of g-mode pulsations studied with single-site data. The noise 
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Table 2. Amplitudes and phases of the leastsquares fit to the merged U, B and V+ light curves of HD 180642. 



Freq 


U& 


u 


B & 


V 


W &y 


& V 


(d-') 


Ampl 


Phase 


Ampl 


Phase 


Ampl 


Phase 




(mmag) 


(rad) 


(mmag) 


(rad) 


(mmag) 


(rad) 


/i = 5.48694 


82.7(7) 


2.884(9) 


48.1(5) 


2.88(1) 


38.1(5) 


2.86(1) 


2/, = 10.97388 


8.2(7) 


1.48(9) 


6.6(5) 


1.52(8) 


5.3(5) 


1.4(1) 


3/i = 16.46082 


4.3(7) 


2.0(2) 


4.3(5) 


2.1(1) 


3.6(5) 


1.9(2) 


/j" = 0.30818 


4.9(7) 


0.9(2) 


3.1(5) 


1.1(2) 


1.6(5) 


0.9(4) 


f^' = 7.36673 


4.5(7) 


5.4(2) 


3.1(5) 


5.4(2) 


2.2(5) 


5.7(3) 


Var. Red. 


96.6 


% 


95.3% 


81.8 


;% 


Residual rms (mmag) 


11.4 


7.9 




12.9 



Errors in unit.s of the la.st digit are given in brackets. The amplitude is expressed in millimag, and the phase is given in radians. The adopted reference epoch for 
<f> = 0.0 corresponds to HJD 2447963.2564. For each combined dataset the total variance reduction of the quintu-periodic model and the root mean square (rms) of 
the residuals is reported. 



level in Fourier space must be computed at each prewhitening 
stage, after checking how the noise level had decreased because 
of prewhitening. 

In the combined Stromgren and ultraviolet datasets, how- 
ever, we found evidence of two additional significant fre- 
quencies. In the Stromgren dataset, we detect 0.2965±ld"' 
and 7.3697 d ', while the ultraviolet data infer the frequencies 
1.30818+1 d"' and 7.36673d"'. Similar frequencies are also 
present in the blue and visual datasets, but are of low amplitude. 
The results from the spectroscopic analysis (see Sect. |3.3l l enable 
us to identify the correct alias peak, and we accept the following 
additonal frequencies from the photometric data of HD 180642: 

=0.30818(5)d-i and/f =7.36673(7) d"'. Figure|2]shows the 
SCARGLE periodograms calculated from the combined Geneva 
U and Stromgren u data (ultraviolet). The residual data (bottom 
part Fig.|2]i still show signs of frequencies at low amplitude, but 
to detect them we would need high-accuracy data of longer time 
span. The U, B, and V datasets, after prewhitening with /i , 2/i, 
and 3/i, folded with and are shown in Fig. [T] The ampli- 
tudes and phases for /i , 2/i , 3/i , and f^' of the leastsquares 
fit to the merged U, B, and W+ light curves of HD 180642 are 
listed in Table|2] 

Within their estimated errors, neither nor f^' is detected in 
the CoRoT light curve. The frequency values found in the space 
photometry are = 0.2 9917(9)d-' and f^' = 7.3586(2)d-i 
(Degroote et al. |20091 l. In addition, after prewhitening 
the main peaks, the frequencies f^- - 0.3204(3) d"' and 
f^^ = 7.3744(3) d"' are significant in the CoRoT periodogram. 
We point out that ~ (f^' + f^^-)/2 and f^' ~ (f^' +f^'-)/2. 
Both sets of three close frequencies cannot be interpreted as ro- 
tational multiplets since it would require a configuration such 
that the central peak is seen in the ground-based data and not 
in the space data. We also exclude the possibility of them be- 
ing independent frequencies coincidentally close to each other 
with a largely varying beat phenomenon because it is unlikely 
to happen twice. Another possiblity is that these close frequen- 
cies are attributed to one and the same mode. This interpreta- 
tion is only plausible if both the error estimates of frequencies 
are largely underestimated and the assumption of constant am- 
plitude and phase during prewhitening of /j^' and f^' from the 
CoRoT light curve produces an artificial additional frequency 
peak close by. The analysis of the very high-precision CoRoT 
light curve indeed showed that the amplitudes, the phases, and 



even the frequencies themselves, f^' and f^' , are not necessarily 
constant in time (Degroote et al. 2009), while constant phases 
and amplitudes for each frequency are traditionally assumed to 
estimate the values of these quantities and their errors as well 
as to prewhiten. It is unclear if this variability in the amplitudes, 
phases, and/or frequencies could be linked to the changing posi- 
tion of the satellite with respect to the Sun in its 5-month orbit or 
modes of truly time-variable behaviour In any case, our ground- 
based photometry is not accurate enough to observe any ampli- 
tude and phase variations with time. Its time base (3043 days) is 
much longer than the time base of the CoRoT observations (156 
days). We consequently interpret the two triple frequency sets 
each being caused by a single oscillation mode, and we keep the 
ground-based frequency values to perform the mode identifica- 
tion from these data. 

2.3. Mode identification 

To identify the modes, we applied the usual method of photomet- 
ric amplitude ratios, using the formalism of Dupret et al. ( |2003l l. 
This technique has already been successfully applied to the pres- 
sure modes in /3 Cephei stars (e.g. Aerts et al. 2006) and to the 
gravity modes in slowly pulsating B stars (De Cat et al. l2007l l. 

We computed the non-adiabatic eigenfunctions and eigen- 
frequencies required for modes with £ < 4- hy using the code 
called MAD (Dupret et al. 2001). Ground-based data are cur- 
rently not precise enough to detect higher degree modes since 
cancellation effects leave a very tight signal that is below the 
detection threshold. The non-adiabatic computations were per- 
formed on two grids of theoretical models for which the adopted 
input physics were different, as described below. All the mod- 
els were computed with the Code Liegeois d' Evolution Stellaire 
(CLES, Scuflaire et al. 2008 ). 

The first grid (grid 1 hereafter) is adopted from De Cat et al. 
(120071 1 and we refer to this paper for details of the input physics. 
The second grid (grid 2 hereafter) was computed using different 
opacity tables and different atmosphere models. For grid 2, OP 
opacity tables (Seaton 120051 1, assuming Asplund et al. (|2005l l 
metal mixtures, were chosen instead of OPAL opacity tables 
(Iglesias & Rogers 1996). Grey atmosphere models were used 
instead of Kurucz atmosphere models. No overshooting was as- 
sumed in grid 1 whereas different values of the overshooting pa- 
rameter (between and 0.5 local pressure scale heights) were 
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Fig. 3. Results of the photometric mode identification for /i (top), (middle) and (bottom). The vertical lines indicate the 
observed amplitude ratios computed from the merged U, B, V+ datasets (see Table|2]i, with their standard error The theoretical 
amplitude ratios are represented for modes of different ^-values. 

Table 3. Logbook of the spectroscopic observations of HD 180642 obtained between June 2005 and July 2007 (HJD 2453536.4- 
2454302.8; AT = 967.5 days). 



instrument 


N 


HJD begin 


HJD end 


AT 


<S/N> 


S/N-range 


T 

^ exp 


resolution 


wavelength range 


FEROS 


224 


2453881 


2454302 


967.5 


186 


[100, 300] 


[300, 1200] 


48 000 


[3570, 9215] 


SOPHIE 


24 


2454265 


2454275 


10.1 


106 


[80, 140] 


[1400, 1500] 


40000 


[3873, 6945] 


Aurelie 


14 


2453536 


2453544 


8.1 


93 


[80, 105] 


900 


25 000 


[4530, 4673] 



For each instrument the number of high-quality spectra N (i.e., spectra with a S/N-ratio > 80 calculated in the region near 4550 A), the HJD, the time-span (in days), 
the average S/N-ratio, the range of S/N-ratios, the typical exposure times (in seconds), the resolution of the spectrograph and the wavelength range (in A) are given. 



considered in grid 2. Finally, grid 2 covers higher masses (up to 
18 Mo) than grid 1 (up to 15 Mq). 

The computation of grid 2 was motivated by the following 
reasons. Since the paper by De Cat et al. ( 120071) . work on the 
modeling of /3 Cephei stars pointed out that OP opacity ta- 
bles can account for the excitation of modes observed in stud- 
ied objects, while OPAL opacity tables fail (Miglio et al. 120071 
Briquet et al. 2007, Dziembowski & Pamyatnykh 2008). In ad- 
dition, closer agreement between theory and observations is also 
achieved if overshooting is present (Briquet et al. 120071 Aerts 
I20081) . 

We obtained the same conclusions regardless of the grid 
considered. In the following, we present our mode identifica- 
tion outcome for grid 2. We used the models whose stellar pa- 
rameters are compatible, within a 3cr error, with those deter- 
mined from spectroscopy in Sect. |3.2| T^ff =24 500+ 1000 K and 
log g =3.45±0.15. 

The outcome for /i is displayed in the top panel of Fig. [3] 
It confirms that the dominant mode is a radial one as already 
ascertained by Aerts ( 120001 ). To identify the ^-values of the two 



observed low-amplitude modes, we then considered models only 
in the 3cr error box that also fit /i as a radial mode. The ampli- 
tude ratios of /j'' and f^' are plotted in the middle and bottom 
panels of Fig. [3] respectively. We conclude that is unambigu- 
ously identified as a g-mode pulsation with i = 3. IVIoreover, f^' 
corresponds to ^ = or 3, with a preference for { = 0. 

It is noteworthy that a low frequency such as 0.30818d"\ 
which corresponds to a high-order g mode, is not yet a common 
feature observed for most of the j0 Cephei stars. However, it is 
certainly not exceptional. For several of the /3 Cephei stars that 
have indeed been the subject of intensive observing campaigns, 
low frequencies have been found, e.g., in 19Mon (Balona et al. 
2002 0.17019 d-'),v Eri (H andler et al. l2004l 0.432 1 8 d" ' ), and 
12 Lac (Handler et al. 120061 0.35529d-i). 

3. Spectroscopy 

3. 1. Data description 

In addition to the photometric datasets, we have at our dis- 
posal 262 useful (S/N > 80) high-resolution spectra gathered 
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with three instruments. HD 180642 was observed during 15 
nights (25 June-4 July; 16-20 July 2007) with FEROS at the 
2.2m ESO/MPI telescope at La Silla and during 12 nights 
(13-24 June 2007) with SOPHIE at the 1.93m telescope at 
OHP in the framework of the CoRoT ground-based Large 
Programme (Uytterhoeven & Poretti 120071 Uytterhoeven et al. 
12008 ). Additionally, 22 Aurelie spectra at the 1.52m telescope 
(OHP) and 1 1 FEROS spectra were already taken in 2005 and 
2006 (14, 17, 19, and 22 June 2005 at OHP; 25-26 May 2006 at 
ESO). An overview and logbook of the spectroscopic observa- 
tions are given in Table |3] 

We carefully calculated the integration times that we selected 
to ensure a S/N ratio of about 200 near the wavelength 4550 A, 
which is the position of a prominent Si III triplet well suited to 
our line-profile variability study (Aerts & De Cat 2003). These 
silicon lines are indeed sufficiently strong without being much 
affected by blending. Moreover, they are dominated by temper- 
ature broadening, such that the intrinsic profile can be modelled 
by a Gaussian. This simplifies the modeling of the line-profile 
variations for mode identification purposes (De Ridder et al. 
12002) . 

We reduced the FEROS spectra using an improved version 
of the standard FEROS pipeline, written in MIDAS, developed 
by Rainer (2003). The main improvements of this pipeline con- 
cern the blaze and flat-field correction of the spectra, by using an 
accurate definition of the blaze function extracted from a well- 
exposed spectrum of a hot star. The SOPHIE spectra were ex- 
tracted and automatically reduced in real-time by a reduction 
package adapted from HARPS. Since there is no pipeline re- 
duction available for the Aurelie spectrograph, we used standard 
reduction procedures with IRAF. After a correction to the he- 
liocentric frame, the spectra were manually normalised using a 
cubic spline fit. 

3.2. Abundance analysis 

The non-local thermodynamic equilibrium (NLTE) abun- 
dances of helium and the dominant metals have been cal- 
culated using the latest versions of the line-formation codes 
DETAIL/SURFACE (Butler & Giddings 1985; Giddings 1981) 
and plane-parallel, fully line-blanketed atmospheric models with 
a solar helium abundance (Kurucz 1993 ). Our analysis is based 
on a mean FEROS spectrum created by coadding 1 1 individual 
exposures obtained in May 2006. All exposures were transferred 
to the laboratory rest frame prior to this operation. Curve-of- 
growth techniques were used to determine the abundances using 
the equivalent widths of a set of unblended lines. The reader 
is referred to Morel et al. (120061 ) for complete details of the 
methodology used to derive the elemental abundances. 

3.2.1 . Atmospheric parameters 

A standard, iterative scheme was used to self-consistently de- 
rive the atmospheric parameters solely by spectroscopic means: 
Teff was determined from the Si ll/lll/iv ionisation balance, log g 
from fitting the collisionally-broadened wings of the Balmer 
lines, and the microturbulent velocity, ^, by requiring that the 
abundances corresponding to the O II features were indepen- 
dent of the line strength. We obtain: reff=24500±1000 K, 
log ^=3. 45 ±0.15 dex and ^=12±3 km s"'. As can be seen in 
FiglH satisfactory fits are obtained for all Balmer lines con- 
sidered. Within the uncertainties, identical values are obtained 
when analysing these data with the unified code FASTWIND 
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Fig. 4. Comparison between the synthetic (dashed line; 
re(f=24500 K and logg=3.45 [cgs]) and observed (solid line) 
Balmer line profiles. The synthetic spectra have been convolved 
with a rotational profile with v=44 km s"'. No attempts have 
been made to fit the other spectral lines. 



(Puis et al. 120051 ) and (semi-)automatic line-profile fitting tech- 
niques (Lefever et al. |2U09j : reff=24 000± 1000 K, log g=3.4+0. 1 
dex, and ^=13+2 km s"'. Support for our rgft estimate is pro- 
vided by the ionisation balance of neon, which suggests an iden- 
tical value: 7^(1=24 OOOilOOO K (Morel & Butler 2008). Only 
a single C III, N III, or S III line could be measured, but good 
agreement with the abundances yielded by the far more numer- 
ous transitions of lower ionisation ions (namely C II, N II, and 
S II) is also found in each case. The fulfillment of ionisation bal- 
ance for all species under investigation indicates that the effec- 
tive temperature of HD 180642 is well constrained. 

Given the large amplitude of the radial mode, we investigated 
the amplitude of the variation in atmospheric parameters over the 
pulsation cycle by analysing separately the average of 20 FEROS 
exposures corresponding to minimum and maximum EWs of the 
Si III fines, and obtained Areff ~ 1 000 K and A log g~0. 1 dex. The 
same result was already obtained for the /3 Cephei star ^' CMa, 
which resembles HD 180642 from the viewpoint of the domi- 
nant mode (Morel et al. 2006). Finally, we infer a total amount 
of line broadening arising from pulsations and rotation, vy=44 
km s by comparing the profiles of a set of isolated O II lines 
with a grid of rotationally-broadened synthetic spectra. The stel- 
lar parameters suggest that HD 180642 is an object close to the 
end of the core-hydrogen burning phase. We refer to a forthcom- 
ing paper on modeling (Thoul et al., in prep.) for an in-depth 
discussion of its evolutionary status. 

3.2.2. Elemental abundances 

The abundances computed for reff=24500 K, logg=3.45 dex, 
and ^=12 km s ' are given in Table IH and compared with the 
standard solar mixture of Grevesse & Sauval ( 119981 ) and val- 
ues derived from time-dependent, three-dimensional hydrody- 
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Table 4. Mean NLTE abundances (on the scale in which 
loge[H]=12) and total 1-cr uncertainties (in brackets), in units 
of the last digits. 





HD 180642 


B stars 


Sun 1-D 


Sun 3-D 


He/H(4) 


0.088(18) 


-0.089 


0.085(1) 


0.085(2) 


log 6(C) (9) 


8.21(10) 


-8.20 


8.52(6) 


8.39(5) 


log6(N) (21) 


8.00(19) 


-7.79 


7.92(6) 


7.78(6) 


log 6(0) (25) 


8.53(14) 


-8.58 


8.83(6) 


8.66(5) 


log6(Ne) (6) 


7.86(16)° 


-8.07 


8.08(6) 


7.84(6) 


loge(Mg)(l) 


7.34(20) 


-7.48 


7.58(5) 


7.53(9) 


log6(Al) (3) 


6.22(15) 


-6.08 


6.47(7) 


6.37(6) 


log6(Si) (7) 


7.19(19) 


-7.20 


7.55(5) 


7.51(4) 


log6(S) (4) 


7.10(34) 


-7.19 


7.33(11) 


7.14(5) 


log6(Fe) (21) 


7.34(21) 


-7.36 


7.50(5) 


7.45(5) 


[N/C] 


-0.21(22) 


—0.41 


-0.60(9) 


-0.61(8) 


[N/0] 


-0.53(24) 


—0.79 


-0.91(9) 


-0.88(8) 


Z 


0.0099(16) 


-0.0108 


0.0172(12) 


0.0124(7) 



The number of spectral lines used is given in brackets in the first column. For 
comparison purposes, we provide typical values found for early B dwarfs in the 
solar neighbourhood (Morel 2009 and references therein; B stai's), the standard 
solar composition of Grevesse & Sauval ( 1998 Sun 1-D), and updated values 
derived from three-dimensional hydrodynamical models (Asplund et al. 2005 , 
Sun 3-D). We define [N/C] and [N/O] as log[e(N)/e(C)] and log[e(N)/e(0)], 
respectively. The metallicity, Z, is given in the last row. To compute this quantity, 
we assumed the abundances of Grevesse & Sauval 1 1998i for the trace elements 
not under study. 

" A negligible difference amounting to 0.01 dex is obtained when using the T^ff 
derived from the Ne ionisation balance (Morel & Butler l2008t . 



namical models (Asplund et al. 120051 1. The quoted uncertainties 
take into account both the line-to-line scatter and the errors aris- 
ing from the uncertainties in the atmospheric parameters. We 
note that the strong, numerous N II features can also be used to 
constrain the microturbulence, and suggest a significantly lower 
value than the Oil lines: ^=7 km s ' (see, e.g.. Trundle et al. 
12004 for a general discussion of this problem). Adopting this mi- 
croturbulent velocity would lead to a slight (A log e < 0.2 dex) 
upward revision of the abundances and metallicity (see Table |5]l 
bringing them into closer agreement with the probably more re- 
liable solar abundances of Asplund et al. (2005). 

The logarithmic ratios of the CNO abundances [N/C] to 
[N/O] point towards a mild nitrogen overabundance, as unex- 
pectedly observed in other slowly-rotating fi Cephei stars (Morel 
et al. 12006). For HD 180642, however, this result cannot be 
firmly established owing to the relatively large error bars and 
the lack of boron data, which could in principle be used to cor- 
roborate (or refute) the existence of such an N excess (see Morel 
et al. |2008l l. The abundances of the other metals are fully com- 
patible with literature values for nearby main-sequence B stars 
(e.g., Daflon & Cu nha 200 4; Gummersbach et al. [19981 Kilian- 
Montenbruck et al. |1994| l. We note that the global metallicity, 
Z =0.0099+0.0016, is very well constrained, because the species 
analysed account for up to 97% of the total metal content. 

3.3. Frequency analysis 

To increase the S/N of the spectra and improve the detection of 
the line-profile variations, we combined several line profiles as 
follows. First, cross-correlated profiles were computed using the 
leastsquares deconvolution (LSD) method (Donati et al. 1 19991 1 



Table 5. Abundance and metallicity differences when adopting 
the microturbulent velocity yielded by the Nil lines (^=7 km 
s '), instead of the adopted value estimated from the Oil fea- 
tures (^=12 km s"'). 





from 12 to 7 km s ' 


AHe/H 


-1-0.025 


A log e(C) 


-1-0.04 




+0.13 


A log 6(0) 


-hO.16 


A log 6(Ne) 


-^0.03 


Alog6(Mg) 


-hO.15 


Alog6(Al) 


-hO.11 


Alog6(Si) 


-hO.10 


Alog6(S) 


-^0.02 


A log 6(Fe) 


-^0.06 


[N/C] 


-^0.09 


[N/O] 


-0.03 


Z 


-^0.0027 



by means of a mask consisting of lines from the VALD database 
(Pisku nov et al. |1995i Ryabchikova et al. 119991 Kupka et al. 
H999] l for Teff = 24 000 K and log^ = 3.50. All elements were 
taken into account, except He and H. This resulted in combining 
about 570 lines for the FEROS and SOPHIE spectra and only 
78 lines for the Aurelie spectra, leading to average S/N ratios of 
600, 400, and 90, respectively. 

Second, to avoid the effects of the different behaviour of the 
individual lines during the pulsation cycle, which was ignored 
in the first method, we computed mean spectra of sufficiently 
strong and unblended lines originating from the same element 
and ionisation state. The averaged spectrum was computed for 
He I lines, N II lines, O II lines, and Si III lines. The use of LSD 
profiles or different elements did not reveal additional periodic- 
ities. In what follows, we describe only our frequency analysis 
of the combination of the two deepest silicon lines of the Si III 
triplet around 4567 A, of an average S/N ratio of 350. 

To perform our frequency analysis, we used the software 
package FAMIASQ (Zima 2008). We first examined the first 
three moments < v' >, < > and < > (see Aerts et al. 1 19921 
for a definition) of the combined silicon lines. Since the line 
profiles move significantly because of pulsation, the integration 
limits for computing the moments were dynamically chosen by 
sigma clipping to avoid the noisy continuum (see Zima 12008] ! . 

Harmonics up to the 5th of the dominant frequency are of 
significant amplitude in < v' >. After removing them, we do find 
four additional significant frequencies in the residuals of < v' >. 
All frequencies retained, except the last one, have an amplitude 
larger than the 4 S/N-level, computed in the residuals, in a 2 d ' 
interval centered on the frequency of interest. The reason for 
retaining the last peak is that it is a one-day-alias of a combi- 
nation frequency. In our frequency solution, we use the value of 
the combination frequency. The detected frequency values and 
their corresponding amplitudes and phases of the leastsquares 
sine fits to < v' > are listed in Table|6] These frequencies were 
also found in the CoRoT light curve (Degroote et al. 2009), albeit 
not with exactly the same values. We akeady discussed a prob- 

' FAMIAS has been developed in the framework of the FP6 European 
Coordination Action HELAS - http://www.helas-eu.org/, 



8 



Briquet et al.: Ground-based observations of the jS Cephei CoRoT main target HD 180642 



•< 



cn 
E 

A 

'"> 

V 



0.37 

0.35 

0.33 

0.31 

0.29 

40 
20 


-20 
-40 
-60 
5 



E 



- 1 



-V:. 



n r 



0.2 



0.4 



0.6 



0.8 



Phase (fi =5.4869 d"') 

Fig. 5. Phase diagram of the equivalent width of the Si III 4552 A 
Une, and of the first three velocity moments, for/i - 5.4869 d"'. 
The dashed lines represent our best model for the moment vari- 
ations (see Sect. l3.4l for more explanations). 



able explanation of these differences in Sect. |2.2| When com- 
paring with our ground-based photometry, we note that the low- 
frequency mode is also found, while the photometric frequency 
at ~ 7.36 d ' is not detected. Moreover, we detect additional pe- 
riodicities in the spectroscopic dataset. A frequency analysis in 
< > and < > did not detect additional periodicities. 

Phase diagrams of < y' >, < >, and < > for the dom- 
inant mode are shown in Fig |5] Several characteristics are no- 
ticeable in < v' >. As already indicated by the photometry, the 
dominant mode does not behave sinusoidally at all and < v' > 
has a peak-to-peak amplitude of about 90 kms"'. This makes 
HD 180642 the /3 Cephei star with the third largest peak-to-peak 
radial-velocity amplitude, BW Vul having ~ 180km s ' (Crowe 
& Gillet ri989l Aerts et al. |l995t M athiase t al. |T998 ) and crSco 
having ~ llOkms ' (Mathias et al. |1991| l. More unusual is the 



Table 6. Frequencies, amplitudes, and phases of the leastsquares 
sine fits to the observed first moment. 



ID 


Frequency 


Amplitude 


Phase 






(km s-') 


(rad) 


ft 




Jo.o(,Z ) 




2/i 


10.9738(1) 


3.0(2) 


5.45(8) 


3/i 


16.4607(1) 


4.7(2) 


3.30(5) 


4/i 


21.9476(1) 


3.6(2) 


5.30(7) 


5/i 


27.4345(2) 


1.7(2) 


0.8(1) 


f2 


8.4079(2) 


1.5(2) 


4.2(1) 




0.3046(2) 


1.1(2) 


2.0(2) 


n 


7.1037(3) 


1.0(2) 


2.7(2) 




13.8948(3) 


0.8(2) 


6.3(3) 



The error estimates in units of the last digit are given in brackets. The adopted 
reference epoch for cji = 0.0 corresponds to HJD 2447963.2564. 



presence of a so-called stillstand, i.e., a plateau at half the de- 
scending branch of the radial velocity curve. In the ground-based 
photometry, this stillstand phenomenon is well observed at cer- 
tain epochs but is not clearly visible in the phased light curves 
of the merged datasets (Fig.[TJ. Such a stillstand is also encoun- 
tered in BW Vul and cr Sco. These common characteristics of the 
three targets suggest a behaviour of the atmospheric pulsation of 
HD 180642 similar to that of BW Vul and crSco. For the latter 
stars, the picture requires the passage of two shocks per pulsa- 
tion period, the stillstand corresponding to a relaxation phase 
(Mathias et al. fT998l l. 

A first indicator of the shock propagation is the evolution in 
the line profile, which exhibits a doubling phenomenon when 
the front shock enters the line forming region, following the 
Schwarzschild mechanism (Schwarzschild |1952| l. The second 
indicator is a phase lag between line forming regions, due to 
the finite shock front velocity, which is present in progressive 
waves. For HD 180642, no line-doubling behaviour is seen in 
the profiles, but it could be hidden by the presence of the low- 
amplitude modes. No phase lag is also observed between varia- 
tions associated with different ions. We therefore conclude that 
the presence of shocks in HD 180642 is not established, nor can 
be rejected. The observation of the star during a few consecutive 
cycles would allow us to prevent the beat phenomenon with the 
non-radial modes and identify the same behaviour as observed 
for BW Vul and cr Sco. 

The behaviour of < > (see middle panel of Fig.|5]) also 
deserves some comments. It clearly cannot be described by a 
double sine, as linear theory predicts for radial pulsators (Aerts 
et al. |1992| l. In addition, the maxima (minima) at phases 0.05 
and 0.45 (at phases 0.25 and 0.80) are different. This difference 
cannot be attributed to a non-linear velocity at the stellar surface 
only but it reflects that temperature variations also take place. By 
taking into account temperature effects, < > can be satisfacto- 
rily modelled, as shown in Sect. 13.41 This is because the second 
moment depends on the thermal width of the local profile, which 
varies with temperature. Similar behaviour, but less prominent, 
was already observed for 6 Ceti (Aerts et al. 1 19921 1 and ^' CMa 
(Saesen et al. l2006] l. 

The temperature variations caused by the large-amplitude 
radial mode are also revealed in the equivalent width of the 
spectral lines. As shown in Fig.|5] the equivalent width of the 
Si III 4552 A line varies with the dominant frequency and two 
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Fig. 6. Amplitude and phase distributions (points with error bars) across the combination of the two deepest silicon lines of the Si III 
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significant harmonics. We point out that this strong EW vari- 
ability is not common among pulsating B-type stars. We refer 
to both De Ridder et al. (2002) and De Cat (2002) for typical 
examples of jS Cephei stars and slowly pulsating B stars, respec- 
tively. However, this behaviour is certainly not unique, and is ob- 
served in several stars with a high-amplitude mode. For /3 Cephei 
stars, it is also the case for BW Vul (Mathias et al.[1998), cr Sco 
(Mat hias e t al. |19911 l, 6 Ceti (Aerts et al. J992] l, v Eridani (Aerts 
et al. l2004l) . and the prototype /3 Cephei itself (Telting et al. |19971 
Catanzaro & Leone I2008D . All these targets are monoperiodic 
radial pulsators or multiperiodic pulsators with a highly domi- 
nant radial mode, as is the case for HD 180642. A different case 
concerns the /3 Cephei star 12 Lac (Desmet et al. ,2009j . This lat- 
ter object, which also is of relatively high amplitude, shows an 
equivalent width that varies with its two dominant modes but the 
highest amplitude mode is identified as an {£, m) = (1,1). 

A frequency search was carried out for the spectra by means 
of a two-dimensional Fourier analysis available in FAMIAS. 
Once the frequencies are selected, a non-linear multi-periodic 
leastsquares fit of a sum of sinusoidals is computed with the 
Levenberg-Marquardt algorithm. This fitting is applied for ev- 



ery bin of the spectrum separately according to the formula 
Z+Yii Ai sin|^27r(/jf -i-(/>,)j, where Z is the zeropoint, and A,-, fi, and 
are the amplitude, frequency, and phase of the /-th frequency, 
respectively. The amplitude and phase distributions across the 
combined silicon lines are shown, in Fig.|6] for the dominant ra- 
dial mode and, in Fig.|2l for the other frequencies detected. 



We emphasize that the 2D frequency analysis detected one 
additional frequency, which was also observed in the CoRoT 
photometry {f^ = 6.3248 d '), compared to the ID analysis of 
the moments, which are integrated quantities. Our Fourier analy- 
sis detected a one-day-alias of and we needed the CoRoT data 
to pinpoint the true periodicity. Additional frequencies were also 
discovered by a 2D frequency search in the case of several other 
P Cephei stars (see Telting et al. 119971 for B Cephei; Schrijvers 
et al. 2004 for e Centauri; Briquet et al. 120051 for 6 Ophiuchi), 
pointing out the necessity of performing this kind of analysis, in 
addition to a frequency search in the radial velocity or higher- 
order moments. 
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3.4. Mode identification 

The moment method (Briquet & Aerts l20()3T l and the Fourier pa- 
rameter fit (FPF) method (Zima 2006) are spectroscopic mode 
identification techniques appropriate to identifying both the de- 
gree { and the azimuthal order m of the oscillations of main- 
sequence pulsators hotter than the Sun. Several successful ap- 
plications, in particular to /3 Cephei stars, are available in the 
literature (e.g. Briquet et al. 12005 1 for 9 Oph, Mazumdar et al. 
12006 for y6 CMa, Desmet et al. |2009| for 12 Lac). 

As explained and illustrated in Zima (2006|, the FPF tech- 
nique, based on mono-mode line profiles, cannot be applied in 
case of a large pulsation velocity relative to the projected rota- 
tional velocity, i.e., if the radial velocity amplitude is above 0.2 
vsin/, as for HD 180642. Instead of using amplitude and phase 
distributions for mono-mode line profiles, we thus considered 
the same quantities but computed them for a huge grid of mul- 
tiperiodic line profile time series. We fixed the dominant mode 



as radial and its velocity field at the stellar surface was modelled 
as the superposition of the velocity fields of five linear radial 
modes. 

To model the local temperature and surface gravity varia- 
tions, we adopted the following empirical approach. At each 
time of observation, the equivalent width (EW) and the width 
cr of an intrinsic Gaussian profile were chosen to match the ob- 
served EW variations and second moment variations. We point 
out that this model assumes that the temperature and gravity vari- 
ations are due solely to the dominant mode, which is a reason- 
able assumption in view of the very low amplitudes of the other 
modes. The theoretical moment variations in this simple model 
are shown in Fig.|5] The amplitude and phase distributions are 
displayed in Fig.|6] Our modeling of the non-linear radial mode 
reproduces the observed behaviour quite well, especially for the 
first four harmonics. 
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Table 7. The bestfit solutions of the spectroscopic mode iden- 
tification for the mode with frequency 8.4079 d ' determined 
by the adapted discriminant E, based on the definition in Aerts 
(fT996] l. 



(A|m|) 


i 


V sin ; 






^t.max 




(3,2) 


81.5 


38.5 


77.35 


30.4 


3.0 


0.93 


(3,1) 


30.5 


39.5 
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The inclination angle i is expressed in degrees; v sin i is the projected rotational 
velocity, expressed in km s"' ; Ap is the amplitude of the radial part of the pul- 
sation velocity, expressed in km s"'; and Vr.max and Vt,max are, respectively, the 
maximum radial and tangential surface velocity due to the inode, expressed in 
km s"' . 



To identify the wavenumbers {[,m) of the low-amplitude 
modes, we added their pulsational velocity field to the one of 
the radial mode modelled as described above, and, using a 
value as in Zima (2006), we searched for the solutions that fit 
the observed amplitude and phase distributions computed from 
multi-mode line profiles best. Unfortunately, no good fit could 
be achieved for any of the modes. Next, we used the moment 
method for multiperiodic stars (version of Briquet & Aerts , 2003 J 
also adapted to adequately model the dominant radial mode, but, 
again, the outcome was inconclusive. 

Therefore, we tried to identify the modes individually, fol- 
lowing the version of the moment method of Aerts ( 1996 ). Since 
the width cr of the intrinsic profile is not constant in time, in 
contrast to the assumption in Aerts (119961 ). we defined a slightly 
different discriminant, in which the amplitudes involving cr are 
not taken into account. We also omitted the constant term of 
the second moment because it contains the contribution of all 
the modes. Therefore, the discriminant used is the one of Aerts 
( I1996I I but only with the amplitudes denoted by C, D, F, and G 
(see Aerts 1996^ being taken into account. 

We varied the free parameter of the projected rotation veloc- 
ity, vsini, from 1 to 40 km s"' with a step 0.1 km s"', and the 
inclination angle of the star, /, from 3 to 90° with a step 0.5°. For 
each combination {{, m, i), the amplitude of the mode Ap (given 
by definition of Briquet & Aerts 2003 ) was not a free parameter 
but was determined to fit the observed amplitude of the first mo- 
ment, within its error. For the low-frequency mode, we fixed £ 
to be 3, as derived from the photometry. Finally, the criterion in 
Sect. 2.5 in Briquet & Aerts (2003) excluded i > 4, except (£, m) 
= (4,0), for the mode with frequency 8.4079 d"'. For the other 
modes, this criterion limited neither { nor m. 

We obtained a safe mode identification for the mode with 
frequency 8.4079 d ' only. The discriminant values for the other 
modes turned out to be too high for us to be confident about their 
outcome. The discriminant values for the identified mode are 
given in Table|71 It does not allow us to determine the sign of m. 
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Fig. 8. Histogram for the equatorial rotational velocity of 
HD 180642 derived from the moment method outcome, using 
the discriminant as a weight. 

However, prograde and retrograde modes can be distinguished 
by means of a blue-to-red and red-to-blue respectively descent 
of the phase in the profile. From Fig.|7] we deduce that the mode 
is prograde and denote it with a positive m-value. From Table|7] 
we conclude that the mode is unambiguously identified as {£, m) 
= (3,2). 

By using the discriminant E as a weight, we constructed his- 
tograms for i, vsinj, and the equatorial rotational velocity Vgq, 
as in Mazumdar et al. ( 120061 ). We computed them by consider- 
ing only solutions with {£, m) = (3,2). For / and v sin /, the dis- 
tributions are almost flat. However, Vgq can be constrained (see 
Fig. [8]). By calculating a weighted mean and standard deviation, 
we obtained Vgq - 38+15 km s"'. 

4. Summary 

The yS Cephei star HD 180642 was observed by the CoRoT 
satellite during a run of 156 days in 2007. The space white- 
light photometry revealed the rich frequency spectrum of the 
star (Degroote et al. |20091 ). In the present study, we have 
provided additional information about the target, based on 
both ground-based multi-colour photometry and high-resolution 
spectroscopy. We placed our object in the (Teff, log g) diagram. 
In addition, we derived the chemical abundances of several el- 
ements as well as the metallicity of HD 180642. Finally, we 
placed constraints on the identification of some modes. 

From the spectroscopic data, we determined Teff = 
24500±1000 K and log^ = 3.45+0.15 dex. In addition, a de- 
tailed NLTE abundance analysis showed that the considered 
abundance values are compatible with values found for B stars in 
the solar neighbourhood, except for a mild nitrogen excess. This 
nitrogen overabundance was also discovered for several other 
yS Cephei stars (Morel et al. i2006l |2008i ). The deduced metallic- 
ity forHD 180642 is Z = 0.0099+0.0016 (or 0.0126+0.0016 de- 
pending on the choice of microturbulence). 

In the photometry, three pulsation frequencies were found to 
be significant: 5.48694 d"', 0.30818 d"', and 7.36673 d"'. The 
first mode is highly dominant with an amplitude in the U-filter 
about 15 times larger than that of the other modes. By means of 
the method of photometric amplitude ratios, we unambiguously 
identified the degree of the first two modes as ^ = and { - 
3. For the frequency 7.36673 d ', two possibilities remained, 
namely ^ = or 3. 
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In the radial velocity measurements, the dominant radial 
mode is highly non-linear, presents a so-called stillstand and has 
a peak-to-peak amplitude of ~90 km s ' . This behaviour is typi- 
cal of the presence of shocks. However, our dataset did not allow 
us to prove (or definitely refute) their existence in HD 180642. 

We attempted mode identification for several low-amplitude 
modes found in the spectroscopy. In our first trials, we used mul- 
tiperiodic versions of the moment method and FPF method. Both 
techniques were adapted and implemented to adequately model 
the dominant radial mode but we failed for the low-amplitude 
modes. However, a successful outcome for one of the modes 
was achieved by means of the discriminant of Aerts (1996), 
which we adapted slightly to our case. This discriminant iden- 
tifies the modes independently. In the presence of a highly dom- 
inant mode, it is thus more robust than the multiperiodic discrim- 
inant of Briquet & Aerts (120031 ). We found the frequency 8.4079 
d"' to correspond to a (£, m) - (3, 2). Finally, we derived a value 
for the equatorial rotational velocity of Veq = 38±15 km s"^ 

All these observational constraints, together with the CoRoT 
results, will be used to compute stellar models of HD 180642. 
This will be presented in a forthcoming paper 
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